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1 .  INTRODUCTION 

Fast-neutron  spectrum  measurements  are  important  to  studies  of  radiation  vulnerability  of  Army 
electronic  systems.  The  Harry  Diamond  Laboratories  (HDD  is  presently  investigating  the  Bennett'  neutron 
spectrometer,  which  determines  the  neutron  spectrum  by  analysis  of  recoil  protons  in  a hydrogen  (H2)  or 
methane  (CH4)  gas.  The  spectrometer  is  controlled  by  a minicomputer. 

In  data  analysis,  corrections  must  be  made  for  carbon  recoils  and  wall-and-end  effects.  The 
magnitudes  of  these  corrections  are  not  usually  large,  but  are  time-consuming  to  make.  It  was,  therefore, 
decided  to  develop  a quick  method  of  obtaining  preliminary  results  without  these  corrections.  In  this 
manner,  the  trend  of  the  data  could  be  followed  as  the  data  were  taken,  and  preliminary  spectrum  results 
could  be  evaluated  before  the  time-consuming  final  analysis  of  the  data  was  undertaken. 


2.  STATEMENT  OF  PROBLEM 

When  neutrons  below  1 4 MeV  scatter  from  a detector  with  a gas  containing  protons,  the  energy 
spectrum  of  recoil  protons  is  related  to  the  energy  spectrum  of  incident  neutrons  by 

E 

<frp(E)  = N * T fom  4>n(E0  dE,  (1) 

0 E' 

where  N is  the  number  density  of  hydrogen  atoms  in  the  detector  gas,  T is  the  collection  time  of  the 
experiment,  E is  energy,  and  o(E)  is  the  neutron-proton  elastic  scattering  cross  section.  The  proton  flux  is 
$p(E)  and  is  the  neutron  flux.  It  is  necessary  to  differentiate  equation  (1 ) to  obtain  the  desired  neutron 
energy  distribution. 

As  stated  in  the  introduction,  the  various  sources  of  background  noise  which  are  corrections  to 
equation  ( 1 ) are  ignored  at  this  stage  of  analysis.  Of  course,  these  corrections  must  eventually  be  made,  but 
they  are  sufficiently  time-consuming  that  the  final  results  are  not  available  for  several  weeks  after  the  data 
have  been  taken. 


3.  METHOD  OF  SOLUTION 

The  proton-recoil  spectrometer  is  controlled  by  an  ND81 2 minicomputer.  Since  the  minicomputer  is 
always  at  the  location  where  the  data  are  taken,  it  is  advantageous  to  use  it  for  preliminary  analysis.  An 
assembly-language  computer  program  to  solve  equation  (1 ) for  4>n(E)  was  written  for  the  ND81 2.  It  has  the 
following  features: 

(a)  Since  the  data  are  acquired  and  analyzed  by  the  same  computer,  there  is  no  need  to  transfer 
data  from  one  machine  to  another.  Since  the  only  output  device  from  the  ND81 2 (compatible  with  other 
computers)  is  paper  tape,  any  transfer  of  data  would  be  time-consuming. 

(b)  Since  the  ND81 2 is  owned  by  HDL,  there  is  no  data-processing  expense. 


' E F Bennett  and  T I Yule,  Techniques  and  Analyses  of  Fast-Reactor  Neutron  Spectroscopy  with  Proton-Recoil  Propor- 
tional Counters,  Argonne  National  Laboratory,  ANL-7763 11971 ). 
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(c)  The  system  provides  a visual  display  of  both  neutron  and  proton  spectra  on  an  oscilloscope.  A 
listing  of  the  data  is  also  available. 

(d)  The  analysis  is  done  in  a matter  of  minutes,  so  that  any  problems  associated  with  the  data  may 
be  seen  as  the  data  are  taken,  rather  than  weeks  later.  If  there  are  problems,  adjustments  may  be  made 
immediately,  without  wasting  days  of  valuable  reactor  running  time. 


4.  PROGRAM  DESCRIPTION 


The  equation  relating  the  neutron  spectrum  <J>N(E)  to  the  derivative  of  the  proton-recoil  distribution 
d<{>p(E)/dE  is 


4 (E)  = _ _L_  _L  d(VE) 

9n'  n*T  O(E)  dE 


(2) 


The  data  4>P(E)  consist  of  the  number  of  recoil  protons  detected  for  each  of  several  energy  intervals  of 
equal  size.  Numerical  differentiation  is  a noise-generating  process,  and  the  data  contain  statistical 
fluctuations,  so  a smoothing  process  is  used  to  obtain  the  best  results.  The  method  used  is  that  of  Bennet,2 
who  suggests  that  the  slope  of  a least-squares  line — fit  through  the  data  surrounding  a data  point — be  used 
as  the  slope  of  the  data.  The  width  of  the  interval  over  which  the  least-squares  line  is  fitted  is 

STHW  = FUDG*  7WI**2  + CON/E  (3) 


where  STHW  is  the  half-width  of  the  slope-taking  interval.  Wl  is  the  intrinsic  resolution  of  the  detector, 
and  the  constant,  CON,  adds  a factor  to  account  for  the  statistical  nature  of  the  gas  amplification.  FUDG 
is  an  arbitrary  scale  factor  which  allows  more  smoothing.  Bennett2  demonstrates  that  a FUDG  factor  as 
large  as  3.0  does  not  significantly  affect  the  area  under  a peak,  although  significant  line  broadening  does 
occur  at  that  value.  A value  of  2.0  to  3.0  for  FUDG  is  found  to  give  the  best  compromise  between 
smoothing  and  loss  of  resolution. 

The  program  (listed  in  app  A)  is  written  in  BASC-1 2,  the  assembly  language  supplied  by  the  computer 
manufacturer.1  The  program  consists  of  a main  control  routine,  which  calls  various  subroutines  as  they  are 
needed,  and  the  subroutines.  A package4  supplied  by  the  manufacturer  is  used  to  do  the  floating-point 
arithmetic  since  hard-wired  floating-point  arithmetic  is  not  available. 

In  general,  the  program  functions  in  the  following  manner.  A subroutine  INPT  is  called  to  obtain  the 
calibration  data  and  constants  needed  to  run  the  program.  The  program  then  locates  the  proton-recoil  data 
in  core  and  a least-squares  line  is  fitted  through  the  intervals  determined  by  equation  (3).  For  each  energy 
interval,  the  slope  is  used  in  equation  (2)  to  calculate  the  neutron  flux  spectrum.  The  results  are  stored  in 
the  minicomputer  memory  and  may  be  viewed  on  the  oscilloscope  or  output  to  Teletype  or  tape  cassette 
by  use  of  the  ND-1 075-01  computer  program.5 


F.  Bennett,  Fast  Neutron  Spectroscopy  by  Proton-Recoil  Proportional  Counters,  Nuclear  Science  Engineering  27 (1967),  16 
1 Nuclear  Data,  Incorporated,  Software  Instruction  Manual,  BASC-1 2 General  Assembler,  Palatine,  IL  (1971). 

4 Nuclear  Data,  Incorporated,  ND8I2  Utilities  Manual,  Palatine,  IL0971),  Ch  9 through  12. 

* Nuclear  Data,  Incorporated,  Software  Instruction  Manual,  ND4420  Single/ Dual  Parameter  Monitor,  Palatine,  IL  (1972) 
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4 . 1 Description  of  Main  Subroutines 


INPT  is  the  routine  which  reads  the  calibration  data,  along  with  N,  T,  Wl,  and  FUDG  (See  eq 
(1 ) and  (2).)  Also,  the  IN  and  OUT  groups  must  be  input  to  INPT.  These  input  numbers  tell  the  computer 
where  in  memory  to  locate  the  proton-recoil  data  and  where  to  put  the  neutron  flux  results. 

INIT  takes  the  calibration  data  and  defines  the  energy  scale. 

IOST  accepts  the  input  and  output  group  numbers  and  translates  them  into  locations  in  core. 

BUFRT  computes  the  slope-taking  half  width  and  transfers  the  appropriate  proton-recoil  data 
to  a buffer  location.  It  is  possible  for  equation  (3)  to  give  a width  which  reaches  beyond  the  data.  If  this 
occurs,  the  slope-taking  interval  is  reduced  so  that  it  reaches  only  to  the  end  of  the  data. 

LSR  is  the  least-squares  fitting  routine.  The  sums  necessary  for  the  fit  are  accumulated  in  integer 
(as  opposed  to  floating-point)  format  to  save  running  time.  A special  routine,  TPCHK,  converts  the  sums 
necessary  for  least-squares  fitting  from  integer  to  floating-point  format,  since  they  may  exceed  the  normal 
double-precision  24-bit  length  of  an  integer  in  the  computer.  The  result  of  program  LSR  is  the  slope  of  the 
least-squares  line. 

OPT  computes  the  flux  4>n(E)  from  the  slope  using  equation  (2).  The  result  is  stored  as  an 
integer  in  the  memory  of  the  computer,  where  it  may  later  be  accessed  by  the  ND-1 075-01 . 

4 . 2 Use  of  Program 

First,  the  ND-1075-01  program  must  be  loaded  into  core.  This  program  is  used  to  load  the  raw 
proton-recoil  data.  The  data  may  be  loaded  into  groups  of  any  size,  as  the  derivative  program 
automatically  adjusts  itself  to  fit  the  group  size  specified  by  the  ND-1 075-01 . However,  the  data  may  not 
be  loaded  into  the  buffer  memory  extension,  since  the  derivative  program  has  access  only  to  fields  02  and 
03  of  core. 


Next,  the  floating-point  package  is  loaded  into  field  01,  along  with  the  SQRT  and  operate 
instruction  overlays. 

The  derivative  program  is  loaded  in  field  01 . The  starting  address  is  4000  F01 . Field  00  is  not 
altered,  so  that  the  functions  supported  by  the  1075-01  located  in  field  00  are  not  affected.  These 
functions  include  all  display-  and  data-acquisition  functions  but  not  I/O  routines  or  any  functions  which 
use  the  Teletype. 

A sample  run  of  the  program  is  included  in  appendix  B.  First,  two  points  and  their 
corresponding  energies  are  input.  They  set  the  energy  scale  of  the  data.  Then,  the  count  time  of  the 
experiment,  the  hydrogen  atom  concentration  of  the  proton-recoil  tube,  and  mechanical  width  of  the  tube 
are  input  for  use  in  equations  (1)  and  (2).  Next  are  given  the  input  group  and  output  group,  which  tell  the 
program  where  to  get  and  to  put  the  processed  data.  The  input  group  may  not  equal  the  output  group,  and 
neither  may  access  the  buffer  memory  extension.  Finally,  the  slope-taking  multiplier  is  input  for  use  in 
equation  (3)  . 


7 


pr  — — - ■ — . n i.  iiliiiiiwi^ — . 

^ 


The  input  data  may  take  various  floating-point  forms.  For  example,  the  number  "one"  may  be 
input  as  1.0,  1 , or  1 0E  + 00.  If,  after  a number  is  typed,  the  space  bar  is  pushed  to  terminate  that 
number,  the  program  will  continue  typing  on  the  same  line.  If  a carriage  return  is  pushed,  the  tvping  will 
proceed  from  the  beginning  of  the  next  line. 

When  the  derivative  program  is  finished,  it  returns  control  automatically  to  the  107  *>-01  The 
results  may  then  be  viewed  on  the  oscilloscope.  If  the  derivative  of  another  group  of  data  is  to  be 
determined,  the  1 075-01  must  be  stopped  and  the  address  register  set  again  to  4000  F01  In  this  manner, 
the  derivative  program  may  be  run  any  number  of  times  in  succession  However,  if  the  data  are  to  be 
output,  it  is  necessary  to  reload  and  use  the  1075-01  program.  A typical  running  time  for  the  program  is 
15  to  20  s for  256  points. 


5.  SUMMARY 

A rapid,  assembly-language  computer  program  has  been  written  to  derive  neutron  spectra  from 
proton-recoil  spectra.  The  program  allows  a preliminary  analysis  of  data  to  be  made  without  having  to  rely 
on  an  external  computer.  The  program  is  an  invaluable  aid  to  on-line  evaluation  of  the  quality  ot  data 
during  acquisition. 
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APPENDIX  A.— The  BASC-12  Computer  Program 


This  program  was  written  for  the  ND812  minicomputer  language,  BASC-12,  supplied  by  the 
manufacturer!  The  program  consists  of  a main  control  routine.  The  subroutines  are  called  as  needed. 
Since  hard-wired  floating-point  arithmetic  is  not  available,  a package  is  supplied  by  the  manufacturer  for 
the  floating-point  arithmetic. 


FS  I P = 
Fci.n= 
FSOKT 
FS  IN  = 
FJMP  = 
FNEG  = 
IFIX  = 
EXJK= 
FLCAT 

xmx= 


1007 
1 004 
=-72101 
1006 
6000 
1 003 
72(05 
13  721 
= 72(06 
1131 


/REFERS  TO  ND-1075  GROUP  SIZE 


[FIELD  01 


*852(0 


/MODIFY  FLOATING  POINT  PACKAGE  OUTPUT  ROUTINE  SO  IT  WORKS 
/ 


852t0 

7413 

TCP 

/OUTPUT  CHARACTER 

852)1 

7414 

TOS 

852(8 

6101 

JMP  .-1 

/WAIT  TILL  READY 

*4000 

/CONTROL  PORTION  OF 

PROGRAM 

/THIS 

y 

PORTION  CALLS 

OTHER  ROUTINES  AS  NEEDED 

2(000 

0640 

TWJPS 

2(001 

4071 

INPT 

/INPUT  DATA 

2(002 

0640 

TVMPS 

4003 

4770 

INI  T 

/PREPARE  EFP 

4004 

0 640 

TWJPS 

4005 

4371 

I OST 

/LOCATE  IN  AND  OUT 

4006 

0504 

TULDJ  FO 

4007 

1 131 

XMAX 

/SET  COUNTER 

Nuclear  Data,  Incorporated,  Software  Instruction  Manual,  8ASC- 1 2 General  Assembler,  Palatine,  IL  (1921 1. 
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40 10 

2301 

SUBL  01  /SKIP  FIRST  AND  l-AST  POINTS 

40  i 1 

54  22 

ST  J CNT1 

4012 

1510 

CLR  J 

4013 

2202 

A DDL  02  /FIRST  POINT  WHICH  IS  DONE 

4014 

5420 

ST  J CNT2 

4015 

501  7 

CONj 

LDJ  CNT2  /TAKE  DERI V ONE  POINT  AT  A T 

401  6 

0540 

TWST  J 

401  7 

4567 

PTNO 

4020 

0640 

TWJFS 

4021 

4575 

BUFRT  /PUT  DATA  INTO  BUFFER 

4022 

0640 

TWJPS 

4023 

5146 

LSR  /LEAST  SCUARF.  FIT 

4024 

0640 

TWJPS 

4025 

4473 

OPT  /PUT  FLUX  IN  OUTPUT  GROUP 

4026 

3406 

ISE  CNT2 

4027 

3004 

DSE  CNT1  /LAST  POINT? 

4030 

61  13 

JMP  CON  /NO 

4031 

0604 

TWJMP  FO 

4032 

0200 

0200  /YES, RETURN  TO  1075 

4033 

0000 

CNT1  , 

0000  /COUNTER 

4034 

0000 

CNT2, 

0000  /POINT  NO 

/OUTPUT  MESSAGE  ROUTINE 

4035 

0000 

OMR* 

0000 

403  6 

5301 

LDJ0  OMR  /GET  START  ADRESS  OF  MESSAGE 

4037 

3502 

IS  X OMR 

4040 

541  1 

ST  J MESSG 

4041 

5210 

OMR  1 , 

LDJA  MSS.SG  /GET  CHARACTER 

4042 

1501 

SNX  J /END  OF  MESSAGE? 

4043 

63  06 

CRETURN  OMR  /YES 

4044 

7413 

TCP  /NO,  TYPE  IT 

4045 

7414 

TOS 

4046 

6101 

.IMP  .-1 

4047 

3402 

ISE  MESSG 

4050 

6107 

JMP  OMR 1 /GET  NEXT  CHAR 

4051 

0000 

MESSG  » 

0000 

4052 

0000 

A DDE  , 

OCOO  /WHERE  TO  PUT  DATA 

4053 

1400 

ADR1  , 

1400 

4054 

0000 

STT, 

0000  /INPUT  THRU  FLOAT  POINT  PACK 

4056 

5102 

LDJ  ADR1  /PREP  TO  ENTER  DPP 

4056 

0540 

TV  ST  J 

4057 

0402 

0402 

4060 

0 640 

TVJPS 

4061 

0400 

0400  /ENTER  FPP 

4062 

1400 

FINPUT 

4063 

571  l 

FSTORA  A DDR 

4064 

1001 

FEX I T 

4065 

51  13 

LDJ  A DDR  /NEXT  DATA  GOES  THREE  WORDS 

4066 

2203 

A DDL  03 

10 


A 


'1067 

'1070 


5515 

6310 


ST J ADDR 

cp.eturn  stt 


/WF.  I T 
/GET 
/ 


E TITLES 

DATA  THRU  FLOATING  POINT  PACKAGE 


AO  71 

0000 

INPT, 

0000 

AO  7 2 

O6A0 

Al  , 

TV  JPS 

AO  73 

A035 

OMR 

/OUTPUT  MESSAGE 

AO  7 A 

A232 

TITLE 

/DERIVATIVE  PEOGP. 

AO  7 5 

7103 

XCT  Al 

AO  7 6 

A21  1 

POINT 

/PO I NT  = 

AO  7 7 

5053 

LDJ  DAT1 

/POINT  TO 

A1  00 

5526 

ST  J ADDR 

A101 

6525 

JPS  STT 

/GET  POINT 

A102 

7110 

XCT  Al 

/PRINT  MESSAGE 

A103 

A22 1 

ENE 

/ENERGY  = 

A10A 

6530 

JPS  STT 

/GET  ENERGY 

A105 

71  13 

XCT  Al 

A106 

A21  1 

POINT 

/POINT* 

A10  7 

6533 

JPS  STT 

AJ  10 

71  1 6 

XCT  Al 

A1  1 1 

A22 1 

ENE 

/ENERGY* 

A1  12 

6536 

JPS  STT 

A1  13 

7121 

XCT  Al 

A11A 

A257 

CTM 

/COUNT  TIME 

A1  15 

65A1 

JPS  STT 

A1  1 6 

71  2A 

XCT  Al 

A1  1 7 

A27A 

H2C0 

/H  CONCENTRATION 

A120 

65  AA 

JPS  STT 

A121 

71  27 

XCT  Al 

A122 

A306 

ME  Vi 

/ME  WIDTH 

A123 

65A7 

JPS  STT 

A1  2 A 

7132 

XCT  Al 

A125 

A322 

DFRM 

/DATA  FROM  GRP 

A126 

65  52 

JPS  STT 

A127 

7135 

XCT  Al 

A1  30 

A335 

DTO 

/DATA  TO  GROUP 

A1  3 1 

6555 

JPS  STT 

A1  32 

71  AO 

XCT  Al 

A1  33 

A35A 

FURG 

/MULTIPLIER 

A1  3 A 

65  60 

JPS  STT 

A1  35 

5162 

LDJ  ADR1 

/PREP  TO  : 

A1  36 

05  AO 

T VST  J 

A1  3 7 

0A02 

0A02 

A1A0 

06  AO 

TV  JPS 

A1  A 1 

OAOO 

OAOO 

A1  A2 

5036 

FI.OAD  GF 

/CHANGE  G: 

A IAS 

7405 

I FIX 

A!  AA 

5A3A 

FSTOR  GF 

Al  A5 

503/' 

FI.OAD  GT 

A1  A6 

7405 

IF  IX 

A1A7 

5A3A 

FSTOR  GT 

START  OF  DATA  BUFFER. 


ENTER  FPP 


GFj  GT  TO  INTEGERS 


11 


APPENDIX  A 


4150  1001  FEXIT 

4151  6360  C RETURN  INPT 


/END  OF  BLOCK  1 
// 

IC  AT  4153 
IC  AT  4153 


41  52 

41  53 

DAT  1 , 

PI 

41  53 

0000 

PI  , 

0000 

41  54 

0000 

0000 

4155 

0000 

0000 

4156 

0000 

El  , 

0000 

4157 

0000 

0000 

41  60 

0000 

0000 

41  61 

0000 

P2  t 

0000 

41  62 

0000 

0000 

41  63 

noon 

0000 

41  64 

0000 

F.2, 

0000 

41  65 

0000 

0000 

41  66 

0000 

0000 

41  67 

0000 

CT, 

0000 

41  70 

0000 

0000 

41  71 

0000 

0000 

41  72 

0000 

HC, 

0000 

41  73 

0000 

0000 

41  ”T4 

0000 

0000 

4»  75 

0000 

MU  t 

0000 

41  76 

0000 

0000 

41  77 

0000 

0000 

4200 

0000 

Gr, 

0000 

4201 

ooco 

0000 

4202 

0000 

OQOO 

4203 

0000 

GT  » 

0000 

4204 

0000 

0000 

4205 

0000 

0000 

4206 

0000 

FDGE  t 

0000 

4207 

0000 

0000 

4210 

0000 

0000 

421  1 

0320 

P0r'1T, 

0320 

4212 

0317 

0317 

4213 

0311 

0311 

4214 

0316 

0316 

4215 

0324 

0324 

4216 

0240 

0240 

421  7 

0275 

0275 

4220 

nooo 

0000 

4221 

0305 

ENE, 

0305 

4222 

0316 

0316 

4223 

0305 

0305 

L. 


/INPUT  DATA  BUFFER 


/COUNT  TIME 
/H  CONCENTEATI ON 
/MECH  WIDTH 
/INPUT  GROUP 
/OUTPUT  GROUP 
/WI  MULTIPLIER 

/POINT 

y! 

/ENERGY 

I 

12 

1 


4224 

0322 

0322 

4225 

0307 

0307 

4226 

0331 

0331 

4227 

0240 

0240 

4230 

0275 

0275 

4231 

0000 

0000 

4232 

0304 

TITLE, 

0304 

4233 

0305 

0305 

4234 

0322 

0322 

4235 

0311 

031  1 

4236 

0326 

0326 

4237 

0301 

0301 

4240 

0324 

0324 

4241 

0311 

031  1 

4242 

0326 

0326 

4243 

0305 

0305 

4244 

0240 

0240 

4245 

0320 

0320 

4246 

0322 

0322 

4247 

0317 

0317 

4250 

0307 

0327 

4251 

0322 

0322 

4252 

0301 

0301 

4253 

0315 

0315 

4254 

0215 

0215 

4255 

0212 

0212 

4256 

0000 

0000 

4257 

0303 

CTM, 

0303 

4260 

0317 

0317 

4261 

0325 

0325 

4262 

0316 

0316 

4263 

0324 

0324 

42  64 

0240 

0240 

4265 

0324 

0324 

4266 

031  1 

031  1 

4267 

0315 

0315 

4270 

0305 

0305 

4271 

0240 

0240 

4272 

0275 

0275 

4273 

0000 

0000 

4274 

0310 

H2C0, 

0310 

42  75 

0240 

0240 

4276 

0303 

0303 

4277 

0317 

031  7 

4300 

0316 

0316 

4301 

0303 

0303 

4302 

0305 

0305 

4303 

0316 

031  6 

4304 

0275 

0275 

4305 

0000 

0000 

L 


/TITLE 


/COUNT  TIME 


/H  CONCEN 


1 3 
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4306 

0315 

MEW* 

0315 

4307 

0305 

0305 

4310 

0303 

0303 

431  1 

0310 

0310 

4312 

0240 

0240 

4313 

0327 

0327 

4314 

0311 

031  1 

4315 

0304 

0304 

4316 

0324 

0324 

4317 

0310 

0310 

4320 

0275 

0275 

4321 

0000 

0000 

4322 

0311 

DFPM, 

0311 

4323 

0316 

0316 

4324 

0240 

0240 

4325 

0307 

0307 

4326 

0322 

0322 

4327 

0317 

0317 

43  3 0 

0325 

0325 

4331 

0320 

0320 

4332 

0240 

0240 

4333 

02  75 

02  75 

4334 

0000 

0000 

4335 

0317  DTO* 

0317 

4336 

0325 

0325 

4337 

0324 

0324 

4340 

0000 

00 

4341 

0325 

0325 

43  42 

0324 

0324 

4343 

024'' 

0240 

4344 

0307 

0307 

4345 

0322 

0322 

434  6 

0317 

0317 

4347 

0325 

0325 

4350 

0320 

0320 

4351 

0240 

0240 

4352 

0275 

0275 

4353 

0000 

0000 

/ME  WIDTH 


/INPUT  GROUP 


/OUT  GROUP 


4354 

0327  F’JDG* 

0327 

4355 

031  1 

031  1 

4356 

0240 

0240 

4357 

0315 

0315 

43  60 

0325 

0325 

4361 

0314 

0314 

43  62 

0324 

0324 

4363 

0275 

0275 

43  64 

0240 

0240 

4365 

0000 

0000 
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/END  OF  BLOCK  2 
// 

IS  C%  AT  4366 

4366  OOOO  ;% 

/TRANS  1 -ATE  IN  GP.OUP  AND  OUT  GROUP  NUMBERS 
/ INTO  ACTUAL  LOCATION  IN  CORE 
/ 

4367  0604  RETN , Tl.'JMP  FO 

4370  0200  0200 

4371  0000  I OST  , 

4372  0500 

4373  4202 

4374  0510 

4375  4205 

4376  1131 

4377  1501 

4400  6111 

4401  0514 

4402  1131 

4403  1 604 

4404  0500 

4405  4202 

4406  2301 

4407  1161 

4410  1000 

4411  1302 

4412  0540 

4413  4572 

4414  1604 

44)5  1604 

44*6  0550 

441 7 4573 

4420  0514 

4421  1131 

4422  1 604 

4423  0500 

4424  4205 

4425  2301 

4426  1161 

4427  1000 

4430  1 302 

4431  0540 

4432  4462 

4433  1 604 

4434  1 604 

4435  0550 

4436  4463 

44  3 7 6 3 4 6 


/CONSTANTS  FOR  HYDROGEN  CROSS  SECTION 
4/;40  0016  CONI,  0016  /11010. 


0000  /LOCATE  GROUPS  IN  CORE 

TVLDJ 

GF+2 

TVLDK  /IN  GROUP  = OUT  GROUP? 

GT+2 
NS  JK  J 
SN7.  J 

JMP  RETN  /YES, RETURN 

TVLDK  FO  /NO, SET  INPUT  DATA 

XMAX 

INC  K /NO  OF  PTS  PER  GROUP 

TVLDJ 

GF+2 

SUBL  01  /NO  OF  GROUPS  - 1 

ROTD  J 01  /DOUBLE  PRECISION 

MPY 

LJKFRS 

TVSTJ 

INST  /LOCATION  OF  INPUT  GROUP 
INC  K 

INC  K /DATA  STARTS  IN  FIELD  2 
TVSTK 

INST  + 1 /FIELD 

TW1.DK  FO  /DO  SAME  FOR  OUTPUT  GROUP 

XMAX 

INC.  K 

TVLDJ 

GT+2 

SUBL  01 

ROTD  J 01 

MPY 

LJKFRS 

TVSTJ 

OUST  /LOCATION  OF  OUTPUT  GROUP 

INC  K 

INC  X 

TVSTK 

OUST  + 1 

C RETURN  I OST 
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4441 

2540 

2540 

44  4 2 

2000 

2000 

4443 

0014 

C 0N2 , 

0014 

/4041 . 

4444 

3744 

3744 

444  5 

4000 

4000 

4446 

0014 

C0N3  , 

0014 

/2387. 

444  7 

2251 

2251 

4450 

4000 

4000 

4451 

0010 

C0N4  , 

0010 

/135.5 

44  5 2 

2074 

2074 

44  5 3 

0000 

0000 

4454 

0000 

STR1  , 

0000 

/GENERAL  STORAGE 

4455 

0000 

0000 

44  5 6 

0000 

0000 

4457 

0000 

STR2, 

0000 

/GF.NEL 

STORAGE 

44  60 

0000 

0000 

4461 

0000 

0000 

4462 

0000 

OUST  , 

0000 

/WHERE 

TO  PUT  OUTPUT  DATA 

44  63 

0000 

0000 

/FIELD 

4464 

0000 

FLUX, 

0000 

4465 

0000 

0000 

44  6 6 

0000 

0000 

44  6 7 

1400 

CONS, 

1400 

4470 

0000 

SIG, 

0000 

/CROSS 

SECTION 

4471 

0000 

0000 

4472 

0000 

0000 

/CONVERT  SI, OPE  TO  FLUX  AMO 
/ PUT  INTO  CORRECT  OUTPUT  LOCATION 


44  73 

0000  OPT, 

0000 

4474 

5105 

LDJ  C0N5  /COMPUTE  FLUX 

4475 

0540 

TWSTJ 

44  7 6 

0402 

0402 

4477 

0640 

TWJPS 

4500 

0400 

0400  /ENTER  FLOAT  POINT  PACE 

4501 

0500 

FTWLD  /COMPUTE  H CROSS  SECTION 

4502 

5034 

E 

4503 

4540 

FADD  COM2 

4504 

5530 

FSTOR  STR1 

4505 

0500 

FTWLD 

4506 

5034 

E 

4507 

4536 

FADD  COM4 

4510 

5531 

FSTOR  STR2 

4511 

5151 

FLOAD  CONI 

4512 

6536 

FDIV  STP.l 

4513 

5537 

FSTOR  STR1 

4514 

51  46 

FLOAD  CON 3 

4515 

6536 

FDIV  STR2 

451  6 

4542 

FADD  STr!  1 

4517 

5527 

FSTOR  SIG  /STORE  IT 
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4520 

0500 

FTWLD  /COMPUTE  FLUX 

4521 

5034 

E /FLUX=-(E/(N*SIG*CT)  ) * 

(SLOPE /DEPP) 

4522 

6532 

FDI V SIG 

4523 

0 700 

FTVJMT 

4524 

5416 

SLOPE 

4525 

0640 

FTV.'DV 

4526 

41  72 

KC 

4527 

0640 

FTV.’DV 

4530 

41  67 

CT 

4531 

0640 

FTV.’DV 

4532 

5037 

DEPP 

4533 

1003 

FNEG 

4534 

1007 

FSIP  /FLUX  < 0 7 

4535 

1004 

FCLR  /YES » FLUX  = 0 

4536 

7405 

I FIX 

4537 

5553 

FSTOR  FLUX 

4540 

1001 

FEX I T 

4541 

5017 

LDJ  PT3  /GET  FIELD 

4542 

4557 

ADJ  OUST+1 

4543 

5402 

ST  J PT1 

4544 

5156 

LDJ  FLUX +2  /GET  FLUX  (LOW 

ORDER) 

4545 

0564 

PT1  , 

TWSTJO  FO  /STORE  IT 

4546 

44  62 

OUST 

4547 

3565 

IS?.  OUST 

4550 

6002 

SKIP 

4551 

6410 

JPS  INCR 

4552 

51  65 

LDJ  FLUX+1  /GET  FLUX  (HIGH 

ORDER) 

4553 

7106 

XCT  PT1  /STORE  IT 

4554 

3572 

IS?  OUST 

4555 

6002 

SKIP 

4556 

6403 

JPS  INCR 

4557 

63  64 

P<T3, 

CRETURN  OPT 

45  60 

05  64 

TVJSTJ0  FO 

45  61 

0000 

INCH, 

0000  /INCREMENT  FIELD 

4562 

0340 

TWIS? 

45  63 

4463 

OUST+1 

4564 

3517 

IS?  PT1 

4565 

6304 

CRETURN  INCR 

4566 

7000 

BUFF  , 

7000 

/3UFFER  LOCATION 

45  67 

0000 

PTNO, 

0000 

/NO  OF  POINT  IN  QUESTION 

45  70 

0000 

PFST  , 

0000 

/FIRST  POINT  OP  SLOPE-TAKING  INTERVAL 

4571 

0000 

PTVTH, 

0000 

/HALF  WIDTH 

45  72 

0000 

INST, 

0000 

/DATA  LOCATION  IN  MEMORY 

45  73 

0000 

0000 

/FIELD 

4574 

0000 

CMTR  , 

0000 

/NO  OF  POINTS  TO  TRANSFER 

/TRANSFER  POINTS  IN 

SLOPE- 

TAKING  INTERVAL  INTO  BUFFER 

45  75 

0000 

BUFRT , 

0000 

45  76 

5107 

LDJ  PTNO 

18 


■ 


^ 

APPENDIX  A 


4664 

3574 

IS?.  PFST 

4665 

6002 

SKIP 

4666 

3512 

IS?  PT4 

4667 

3173 

DS?  CNTR 

/LAST  PT? 

4670 

6114 

JMP  PT4 

/MO 

4671 

63  74 

[RETURN  BUFRT 

/YES 

4672 

4661 

ADR5  , 

PT5 

4673 

0000 

PFSET , 

0000  /CONVERT  POINT  NO  TO  CORE  LOCATION 

4674 

1710 

CLR  JK 

4675 

0500 

TVLDJ 

4676 

4570 

PFST  /GET  NUM  OF  FIRST  POINT 

4677 

2301 

SUBL  01 

4700 

1361 

ROTD  JK  01  /DOUBLE  PRECISION 

4701 

1450 

CLR  0 

4702 

0440 

TWADJ 

4703 

4572 

INST  /CONVERT  TO  CORE  LOCATION 

4704 

0540 

TVJSTJ 

4705 

4570 

PFST  /FIRST  CORE  LOCATION  TO  GET 

4706 

1445 

SI?  0 /CHECK  FIELD 

4707 

7315 

XCT6  ADR5 

4710 

1605 

SI?  K 

471  1 

7317 

XCT0  ADR5 

4712 

6317 

C RETURN  PFSET 

/END  OF  BLOCK  3 

// 

IC  AT 

47  i 4 

IS  C 

AT  4713 

4713 

0000 

#; 

4714 

7776 

CONST 1 , 

7776  /0.17 

4715 

2560 

2560 

4716 

5076 

5076 

4717 

0000 

WI  * 

0000  /INTRINSIC  WIDTH  (KEV) 

4720 

0000 

0000 

4721 

0000 

0000 

4722 

0000 

VHP, 

0000  /V.'I  IN  POINTS  (HALF  WIDTH) 

4723 

0000 

0000 

4724 

0000 

0000 

4725 

ooco 

VMSQ# 

0000  /WM*WM 

4726 

0000 

0000 

4727 

0000 

0000 

4730 

1400 

FPC0N1 , 

1400 

4731 

0002 

two. 

0002  /FLOATING  POINT  2.0 

4732 

2000 

2000 

4733 

0000 

0000 

/FIND 

/ 

SLOPE 

TAKING 

HALF  INTERVAL 

4734 

0000 

VIN, 

0000  /VI =SCRT ( WN**2+C0NST1 /E  > 

4735 

5105 

LDJ  FPC0N1 
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4736 

0540 

TWSTJ 

4737 

0402 

0402 

4740 

0640 

TWJPS 

4741 

0400 

0400 

/ENTER 

FLOAT  POINT  PACK 

4742 

5126 

FLO  AD 

C0NST1 

4743 

0640 

FTVDV 

4744 

5034 

E 

4745 

0440 

FT  WAD 

4746 

4725 

WMSO 

474  7 

7401 

FSQPT 

4750 

0700 

FTWMT 

4751 

5034 

E 

4752 

5533 

FSTOR 

WI 

/WI  IN  KEV 

4753 

0640 

FTWDV 

4754 

5037 

DEPP 

/WI  IN 

POINTS 

4755 

0700 

FTWMT 

475  6 

4206 

FDGE 

4757 

6526 

FDI V TWO 

/HALF  INTERVAL 

4760 

7405 

IFIX 

4761 

5537 

FSTOR 

WIP 

/WIDTH  IN  POINTS 

4762 

1001 

FEXIT 

4763 

6327 

ERETURN  WIN 

4764 

0000 

STOP.  1 , 

0000 

4765 

0000 

0000 

/FLOAT 

NO  TEMP  STORAGE 

4766 

0000 

0000 

4767 

1400 

ADR6, 

1400 

4770 

0000 

INIT, 

0000 

/INITIALIZE  SO 

4771 

5102 

LDJ  ADR6 

/FLOATING  CONSTANTS 

4772 

0540 

EN1  , 

TWSTJ 

4773 

0402 

0402 

4774 

0640 

EN2, 

TWJPS 

4775 

0400 

0400 

4776 

0500 

FTWLD 

/FIND 

DEPP  FROM  CALIBRATION 

4777 

41  61 

P2 

5000 

0400 

FTWSB 

5001 

4153 

PI 

5002 

5516 

FSTOR 

ST0R1 

5003 

0500 

FTWLD 

5004 

4164 

E2 

5005 

0400 

FTWSB 

5006 

4156 

El 

5007 

6523 

FDI V ST0R1 

5010 

5427 

FSTOR 

DEPP 

SO!  1 

0500 

FTWLD 

/FIND 

B FOR  E=PT*DEFP+B 

5012 

4153 

PI 

5013 

7024 

FMULT 

DEPP 

5014 

5530 

FSTOR 

ST0R1 

5015 

0500 

FTWLD 

501  6 

41  56 

El 

POINTS 
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501  7 

4133 

FSUB  ST0R1 

5020 

5422 

FSTOR  B 

5021 

0500 

FTV’LD 

/INITIALIZE  WMSQ 

5022 

41  75 

MW 

5023 

0700 

FTWMT 

5024 

4175 

MV? 

5025 

0540 

FTV.'ST 

5026 

4725 

msQ 

5027 

1001 

FEXIT 

5030 

6340 

[RETURN 

INIT 

5031 

0000 

Pi 

0000 

/POINTER  STORAGE 

5032 

0000 

0000 

5033 

0000 

0000 

50  34 

0000 

E* 

0000 

/ENERGY  RESULT 

5035 

0000 

0000 

5036 

0000 

0000 

5037 

0000 

DEPP  * 

0000 

/E  =DEPP  X P + B 

5040 

0000 

0000 

5041 

0000 

0000 

5042 

0000 

B» 

0000 

5043 

0000 

0000 

5044 

0000 

0000 

5045 

0000 

EFP * 0000  /FIND  ? 

ENERGY  FROM  POINT 

504  6 

5157 

LDJ  ADR 6 

5047 

7155 

XCT  EMI 

5050 

7154 

XCT  EM 2 /ENTER 

FPP 

5051 

5120 

FLOAD  P 

5052 

7406 

FLOAT 

5053 

71  14 

FMULT  DEPP 

5054 

4512 

FADD  D 

5055 

5521 

FSTOR  E /E=DF.PP*P+B 

505  6 

1001 

FEXIT 

5057 

6312 

[RETURN  EFP 

/XCLR 

STORES  ZEROS  IN  SUCCESSIVE  LOCATIONS 

/CALL 

s UPS 

XCLR 

/ 

START-  FIRST  V.’OR.D  CLEARED 

/ 

NO 

- NO  OF  V?0RDS  CLEARED 

/ 

J 

RETURN  POINT 

r 

5060 

0000 

XCLR*  0000 

5061 

5301 

LDJ?  XCLR 

5062 

3502 

ISZ  XCLR 

5063 

540  6 

ST J ADI  /FIRST 

UD  TO  BE  ZEROED 

50  64 

5304 

LD JO  XCLR 

5065 

541  1 

ST  J XC2  /NO  OF 

V.'DS  TO  BE  ZEROED 

50  6 6 

3506 

ISZ  XCLR 

/SET  RETURN  ADRESS 

506  7 

1510 

CLR  J 

5070 

0540 

TWSTJ 

5071 

0000 

ADI*  0000 
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5072 

3501 

IS7.  ADI 

50  73 

3003 

DS7.  XC2 

50  7 A 

61  04 

JMP  .-4 

5075 

6315 

[RETURN  XCLH 

5076 

0000  XC2# 

0000 

/END  OF  BLOCK  4 
// 

IC  AT  5100 
5077  0003  #3 

/CONVERT  DOUBLE  OR  TPLE  PRECISION  INTEGER 
/TO  FLOATING  POINT 
/ 


5100 

0000 

TPCHK, 

0000 

5101 

0500 

TWLDJ 

5102 

5407 

POSST 

5103 

5437 

ST J POS 

/LOCATION  OF  NUMBER 

5104 

2302 

SUBL  02 

5105 

5434 

ST  J P0S2 

5106 

0500 

TVLDJ 

/PREPARE  FPP 

5107 

5254 

I DC ONI 

51  10 

0540 

T VST  J 

51  1 1 

0402 

0402 

51  12 

0640 

TWJPS 

5113 

0400 

0400 

/ENTER  IT 

51  14 

0520 

FTVLD& 

/GET  LOW  ORDER  TWO  WORDS  OF  NUMBER 

51  15 

5142 

POS 

5116 

7406 

FLOAT 

51  17 

5424 

FSTOR  ST0R5 

5120 

1006 

FS  IN 

/IF  NEG,  NO  VANTER  =PREC  +N 0 FOUND 

5121 

6004 

FJMP  .+4 

/SINCE  ND  ASSUMES  DOUBLE 

PREC 

INTEGB 

5122 

5014 

FLOAD  PREC  /WITH  HIGH  ORDER  BIT  SET 

TO  BE 

NEG 

5123 

4420 

FADD  ST0R5 

5124 

5417 

FSTOR  ST0R5 

5125 

0520 

FTVLDO 

5126 

5141 

P0S2 

/GET  HIGH  ORDER  PART  OF  MUM 

5127 

7406 

FLOAT 

5130 

7006 

FMULT  PREC  /CORRECT  MAGNITUDE 

5131 

4412 

FADD  ST0P.5  /ADD  LOW  ORDER  PART 

5132 

0560 

FT  VS  T9 

/STORE  FLOATING  NUMBER 

5133 

5142 

POS 

5134 

1001 

FEXIT 

5135 

6335 

[RETURN 

TPCHK 

5136 

0031 

PREC  » 

0031 

/4096*4096 

51  37 

2000 

2000 

5140 

0000 

0000 

5141 

0000 

P052, 

0000 

5142 

0000 

POS, 

0000 

5143 

0000 

ST0R5  p 

0000 

5144 

0000 

0000 

5145 

0000 

0000 
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/LINEAR  LEAST  SOUAEES  FIT  ROUTINE 


/ACCUMULATE 

/ 

SUMS  IN 

INTEGER  FORMAT 

TO  SAVE  TIMS 

/ 

5146 

0000 

LSR, 

0000 

5147 

0640 

TWJPS 

5150 

5060 

XCLR  /CLEAR 

SUMS 

5151 

5231 

I CONI 

5152 

0021 

0021 

5153 

1504 

INC  J / IN I T POINTERS 

5154 

5477 

STJ  ACNT2 

/AND  COUNTERS 

51  55 

5053 

LDJ  NPTS 

5156 

5474 

ST  J ACMT1 

5157 

5472 

ST  J N + 2 

51  60 

0500 

TVLDJ 

51  61 

4566 

BUFF 

51  62 

5400 

ST  J FIRST 

5163 

5070 

C0NTN2, 

LDJ  ACNT2 

/SUM  OF  I 

5164 

1450 

CLR  0 

5165 

4446 

ADJ  ICON  1 +2 

5166 

5445 

ST J ICON  1 +2 

5167 

1445 

SIT  0 /OVERFLOW? 

5170 

3442 

1ST  ICON  1+1 

/YES 

5171 

5062 

LDJ  ACNT2 

/NO*  SUM  I X 

51  72 

1204 

LKFJ 

5173 

1000 

MPY  /I  X I 

55  74 

1302 

LJKFRS 

51  75 

1450 

CLR  0 

5176 

4440 

ADJ  IC0N2+2 

/LOW  ORDER 

51  77 

5437 

S3 J I COM 2 +2 

5200 

1445 

SIZ  0 /OVERFLOW? 

5201 

1 604 

INC  K /YES 

5202 

13  74 

EX JK  /MO 

5203 

4432 

ADJ  I C ON 2 + 1 

5204 

5431 

ST  J IC0M2  + 1 

5205 

5200 

LDJP  FIRST 

/Y  LOW 

5206 

1374 

EX  JK 

5207 

5200 

LD JO  FIRST 

/Y  HIGH 

5210 

1374 

EX  JK 

521  1 

1301 

LRSFJK  /STORE 

DATA 

5212 

1450 

CLR  0 

5213 

4433 

ADJ  Y+2 

/SUM  ON  Y 

5214 

5432 

ST  J Y+2 

5315  1445  SIT.  0 /CHECK  FOR  CARRIES  TO  HIGHER  PRECISION 

5316  3427  1ST  Y+l 

5217  6002  SKIP 

5220  3424  1ST  Y 

5221  1450  CLR  0 

5222  1374  EXJK 

5223  4422  ADJ  Y+l 

5224  5421  STJ  Y+l 

5225  1445  SIT  0 

5226  3416  1ST  Y 

5227  6030  JMP  C0MTM1 
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5230 

0000 

NPTS  > 

0000 

/NO  OF  POINTS 

5231 

0000 

ICON  1 t 

0000 

/SUM  OF  I 

5232 

0000 

0000 

5233 

0000 

0000 

523  A 

0000 

I C ON  2 » 

0000 

/SUM  OF  I X I 

5235 

0000 

0000 

5236 

0000 

0000 

5237 

0000 

0000 

5240 

0000 

IY, 

0000 

/SUM  OF  I X Y 

5241 

0000 

0000 

5242 

0000 

0000 

5243 

0000 

0000 

5244 

0000 

Y, 

0000 

/SUM  OF  Y 

5245 

0000 

0000 

5246 

0000 

0000 

5247 

0000 

N, 

0000 

/NO  OF  POINTS 

5250 

0000 

0000 

5251 

0000 

0000 

5252 

0000 

ACNT1  , 

0000 

/COUNTER 

5253 

0000 

A0NT2, 

onoo 

/I 

5254 

1400 

tDCONl , 

1400 

5255 

0600 

C JMP, 

TWJMP 

5256 

51  63 

C0NTN2 

5257 

1302 

C0NTN1  , 

LJKFRS 

/RESTORE  Y 

52  60 

0550 

TWSTK 

/FIND  I X Y 

5261 

5410 

ST0R6 

/HIGH  ORDER  L< 

5262 

0510 

TWLDK 

5263 

5253 

ACNT2 

/r 

5264 

1000 

MPY 

/I  X LOW  0RDE1 

5265 

1 302 

LJKFRS 

5266 

1450 

CLR  0 

ADJ  IY+2  /ADD  IT  (LOW) 

ST J I Y+2  /STORE  IT 

EX  JK 

SI 7.  0 /CHECK  FOR  CARRIES 

I SZ  IY+1 

SKIP 

ISZ  IY 

CLR  0 

ADJ  IY+1  /ADD  IT  (HIGH) 

ST  J IY  + 1 
5IZ  0 
IS  7.  IY 
TWLDK 
ST0R6 
LDJ  A*CNT2 

MPY  /I  X HIGH  ORDER  - DOUBLE  PRECISION! 

LJKFRS 

CLR  0 

ADJ  IY+1  /ADD  IT  (LOW) 

ST  J IY  + 1 
EX  JK 


/ADD  IT  (LOW) 
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5314 

1445 

SIZ  0 /CHECK  FOR  CARRY 

5315 

1504 

INC  J 

531  6 

4556 

ADJ  IY  /ADD 

IT  (HIGH! 

5317 

5557 

ST  J IY 

5320 

3545 

ISZ  ACNT2 

5321 

3147 

DSZ  ACNT1 

/LAST? 

5322 

7145 

XCT  CJMP 

/NO 

5323 

0640 

TWJPS  /YES.. 

FLOAT  SUM  OF  Y 

5324 

5100 

TPCHK  /CHECK  FOR 

5325 

5062 

LDJ  POSST 

/ TRIPLE  PREC 

5326 

2304 

SUBL  04 

5327 

5460 

ST  J POSST 

/FLOAT  SUM  OF 

5330 

0640 

TWJPS  /MAY 

BE  TRIPLE  PREC I S I 

5331 

5100 

TPCHK 

5332 

5055 

LDJ  POSST 

5333 

2204 

ADDL  04 

5334 

5453 

ST J POSST 

5335 

51  61 

LDJ  IDC0N1 

/START  PART  OF 

5336 

0540 

TWSTJ 

5337 

0402 

0402 

5340 

0640 

TWJPS 

5341 

0400 

0400  /ENTER  FPP 

5342 

0500 

FT  WLD 

5343 

5231 

ICON  1 

5344 

7406 

FLOAT 

5345 

0540 

FTWST 

5346 

5231 

I CONI 

534  7 

0700 

FT  WMT 

5350 

5244 

Y 

5351 

5437 

FSTOR  ST OR 6 

/ST0R6=  ( SUM  IP 

5352 

0500 

FT  WLD 

5353 

5247 

N 

5354 

7406 

FLOAT 

5355 

0540 

FTWST 

5356 

524  7 

N 

5357 

0700 

FT  WMT 

5360 

5240 

IY 

5361 

4027 

FSUB  ST0R6 

53  62 

5431 

FSTOR  ST0R2 

/ST0R2  =N* ( SUM 

5363 

0500 

FT  WLD 

5364 

5231 

I C ON  1 

5365 

0 700 

FT  WMT 

5366 

5231 

I CONI 

5367 

5421 

FSTOR  STOP. 6 

/ST0R6=( SUM  I) 

53  70 

0500 

FT  WLD 

5371 

5234 

I CON  2 

5372 

7406 

FLOAT 

5373 

0700 

FT  WMT 

53  74 

5247 

N 

/ACCUMULATORS 

53  75 

4013 

FSUB  ST0R6 

53  76 

0700 

FT  WMT 

53  77 

5037 

DEPP 

5400 

5410 

FSTOR  ST OR 6 

/ST0R6=DEPP* ( 

I X Y 


ROUTINE  TO  FIND  SLOPE 


I)* (SUM  Y) 


IY)  - ST OR 6 


1*1 >*N 


(sum  i*n*N-(Sun  n**s 
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5401 

5012 

FLOOD  ST0P.2 

540  2 

6406 

FD I V ST0P6 

54  03 

5413 

FSTOR  SLOPE 

5404 

1001 

FEX I T 

5405 

0620 

5406 

5146 

CRETURM  LSR 

5407 

5244 

POSST, 

Y 

5410 

0000 

ST0P6, 

0000  /TEMP 

5411 

0000 

0000 

5412 

0000 

0000 

5413 

0000 

ST0P.2, 

0000 

5414 

0000 

0000 

5415 

0000 

0000 

5416 

0000 

SLOPE, 

0000 

541  7 

0000 

0000 

5420 

0000 

0000 

/SLOPE  = ST  0R2 /ST  OP.  6 


STORAGE 


SE  6203 

A1 

= 4072 

ACNT1 

= 5252 

AC  NT  2 

= 5253 

ADI 

= 5071 

ADDR 

= 4052 

ADR1 

= 4053 

ADR  5 

= 4672 

ADR  6 

= 4767 

B 

= 5042 

BUFF 

= 4566 

BUFRT 

= 4575 

CJMP 

= 5255 

CNT1 

= 4033 

CNT2 

= 4034 

CNTR 

= 4574 

COM 

= 4015 

CONI 

= 4440 

COM2 

= 4443 

C0N3 

= 4446 

COM4 

= 4451 

C0N5 

= 4467 

CONST  1 

= 4714 

C0NTN1 

= 5257 

C0NTN2 

= 5163 

CT 

= 4167 

CTM 

= 4257 

DAT  1 

= 4152 

DEPP 

= 5037 

DFRM 

= 4322 

DTO 

= 4335 

E 

= 5034 

El 

= 4156 

E2 

= 4 1 64 

EFP 

= 5045 
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EN1 

a 

4772 

EN2 

= 

4774 

ENE 

= 

4221 

EX  JK 

3 

1374 

FCLR 

a 

1004 

EDGE 

a 

4206 

FJMP 

a 

6000 

FLOAT 

a 

7406 

FLUX 

a 

4464 

FNEG 

a 

1003 

FPCON1 

s 

4730 

FSIN 

a 

1006 

FS  IP 

a 

1007 

FSQRT 

a 

7401 

FUDG 

a 

4354 

GF 

a 

4200 

GT 

a 

4203 

H2C0 

a 

4274 

UC 

a 

41  72 

ICON1 

a 

5231 

ICONS 

a 

5234 

IDCON1 

a 

5254 

IF  IX 

a 

7405 

INCR 

a 

4561 

INIT 

a 

4770 

INPT 

a 

4071 

INST 

a 

4572 

IOST 

a 

4371 

IY 

a 

5240 

LSR 

a 

5146 

MASK 

a 

4646 

PESSG 

a 

4051 

FEW 

a 

4306 

MW 

a 

41  75 

N 

a 

5247 

NPTS 

a 

5230 

OMR 

a 

4035 

0MR1 

a 

4041 

OPT 

a 

4473 

OUST 

a 

4462 

P 

a 

5031 

PI 

a 

4153 

P2 

a 

4161 

PFSET 

a 

4673 

PFST 

a 

4570 

POINT 

a 

421 1 

POS 

a 

5142 

P0S2 

a 

5141 

POS  ST 

a 

5407 

PREC 

a 

5136 

PT1 

a 

4545 

PT2 

a 

461  7 

PT3 

a 

4560 

PT4 

a 

4654 

PT5 

a 

4661 

PTNO 

a 

4567 

PTVTH 

a 

4571 

HETN 

a 

4367 

SIG 

a 

4470 

SLOPE 

a 

541  6 

ST0R1 

a 

4764 

ST0R2 

a 

5413 

ST  OR  5 

a 

5143 

ST  OR  6 

a 

5410 

STP.l 

a 

4454 

STR2 

a 

4457 

STT 

a 

4054 

TITLE 

a 

4232 

TPCHK 

a 

5100 

TWO 

a 

4731 

VI 

a 

471  7 

VI N 

a 

4734 

VIP 

a 

4722 

VMSQ 

a 

4725 

XC2 

a 

5076 

XCLR 

a 

5060 

XMAX 

a 

1131 

Y 

a 

5244 

ER  0002 
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APPENDIX  B.— SAMPLE  PROGRAM  RUN 


This  appendix  provides  a sample  run  of  the  computer  program  BASC- 1 2. 


DERIVATIVE  PROGRAM 
POINT  =1  ENERGY  =1.1 
POINT  =2.0  ENERGY  =.2E+0t 
COUNT  TIME  = 1 .E+03  H C0NCEN=. ! E-02 
VECH  V I DTH  = * 1 

IN  COUP  =1  OUTPUT  GROUP  =2 
VI  MULT  = 3. 

* DERIVATIVE  PROGRAM 
POINT  =1  ENERGY  =1 
POINT  =2  ENERGY  =2 
COUNT  TIME  =1  H CONCEN= 1 
MECH  V IDTH=* 1 

IN  COUP  =1  OUTPUT  GROUP  =2 
VI  MULT  = 3. 

* 0 
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